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Abstract 

We investigate the motion of neutral test particles in the gravitational field of a mass M with 
charge Q described by the Reissner-Nordstrom (RN) spacetime. We focus on the study of circular 
stable and unstable orbits around configurations describing either black holes or naked singulari- 
ties. We show that at the classical radius, defined as Q'^/M, there exist orbits with zero angular 
momentum due to the presence of repulsive gravity. The analysis of the stability of circular orbits 
indicates that black holes are characterized by a continuous region of stability. In the case of 
naked singularities, the region of stability can split into two non-connected regions inside which 
test particles move along stable circular orbits. 
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I. INTRODUCTION 



In general relativity, the gravitational field of a static, spherically symmetric, charged 
body with mass M and charge Q is described by the Reissner-Nordstrom (RN) metric 
which in standard spherical coordinates can be expressed as 

where A = — 2Mr + Q^, and the associated electromagnetic potential and field are 

A = jdt, F = dA = -^dt Adr , (2) 

respectively. The horizons are situated at r± = M ± y/M'^ — Q'^. 

The study of the motion of test particles in this gravitational field is simplified by the 
fact that any plane through the center of the spherically symmetric gravitational source is a 
geodesic plane. Indeed, it can easily be seen that if the initial position and the tangent vector 
of a geodesic lie on a plane that contains the center of the body, then the entire geodesic 
must lie on this plane. Without loss of generality we may therefore restrict ourselves to the 
study of equatorial geodesies with 9 = 7r/2. 

The tangent vector to a curve a;"(r) is = dx'^/dr = i", where r is an afiine 
parameter along the curve. The momentum = fix"' of a particle with mass fi can be 
normalized so that gapx°'x^ = —k, where k = 0,1, —1 for null, timelike, and spacelike 
curves, respectively. For the RN metric we obtain 

- + ^r^ + r'<p' = -k (3) 

on the equatorial plane. The last equation reduces to a first-order differential equation 

^2^2 ^2 ^2 



yU^A A /i^r^ 

where we have used the expressions for the energy, E = —gapitP^ — A*^^; ^^id angular 
momentum, L = gajs^^P^ = of the test particle which are constants of motion associ- 
ated with the Killing vector fields = dt and ^(/, = d^, respectively. Equation (j4]) can be 
rewritten as 



772 

r^ + V^ = — , with F = , , 



fc + ^Vl-^ + ^V (5) 



The investigation of the motion of test particles in the gravitational field of the RN metric is 
thus reduced to the study of motion in the effective potential V. In this work, we will focus 
on the study of circular orbits for which r = and V = E/fi, with the condition dV/dr = 0. 
A straightforward calculation shows that this condition leads to 

an expression which we substitute in Eq. to obtain (an alternative analysis using an 

orthonormal frame is presented in Appendix 

^ (r2-2Mr + Q2)2 
^2 r2(r2-3Mr + 2Q2) " ^ 

Moreover, from the physical viewpoint it is important to find the minimum radius for stable 
circular orbits which is determined by the inflection points of the effective potential function, 
i.e., by the condition d'^V/d'^r = 0. It is easy to show that for the potential ([5]), the last 
condition is equivalent to 

Mr^ - 6MV2 + 9MgV - AQ^ = . (8) 

In this work, we present a detailed analysis of the circular motion of test particles governed 
by the above equations. We will see that the behavior of test particles strongly depends on 
the ratio Q/M and, therefore, we consider separately the case of black holes, extreme black 
holes and naked singularities. 

II. BLACK HOLES 

From the expressions for the energy and angular momentum of a timelike particle [k = 1) 
we see that motion is possible only for r > Q'^/M = and for — 3Mr + 2Q^ > 0, i. 
e., r < and r > r^^, with r^^ = [3M ± ^ (QM^ — 8Q'^)]/2. In fact, from Eqs.(l6]) and 
([7]) it follows that the motion inside the regions r < and r G (r^_,r^^) is possible only 
along spacelike geodesies. At r = r^^ one finds instead that the velocity of test particles, as 
defined in AppendijXl is Ug^ = 1, i.e., the circle r = r^_^ represents a null hypersurface. On 
the other hand, for a non-vanishing charge in the black hole region one can show that 

r_ < < r* < r+ < r^^, with r+ = r_ = = r* for Q = M (9) 
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FIG. 1: In this graphic the radii r^^ = [3M + V(9M^ - 8Q^)]/2 and r± = M + ^/M'^ - Q'^, and 



Q'^/M are plotted. Timehke circular orbits exist only for r > r-y^, whereas r 



'7+ 



represents 



a null hypersurface. Circular motion inside the regions r < and r G {r-y_ , r-y_^_) is possible only 
along spacelike geodesies. 



and r_ 



for Q = 



(10) 



The location of these radii for different values of the mass-to-charge ratio is depicted in 
Fig. [H In the black hole region we limit ourselves to the study of circular timelike orbits, 
i.e., orbits with r > r^_^. 

The effective potential for a test particle with a fixed value of Q/M is plotted in 
Fig. [21 for different values of the angular momentum L/{M^). At infinity, the effective 
potential tends to a constant which is independent of the value of the parameters of the 
test particle and of the gravitational source. In our case, this constant is normalized by 
choosing the value of the total energy of the particle a.s E/fi. Moreover, as the outer horizon 
is approached from outside, the effective potential reaches its global minimum value which 
is zero. This behavior is illustrated in Fig. |3] where the effective potential is depicted for a 
specific value of Q/M and different values of L/{Mfi). The radius of a circular orbit, rco, 
is determined by the real positive root of the equation 

2Q^L^ 



Mr-" 



Q + — I ^ + —r 



. 



(11) 



In general, in the region r > circular orbits do not always exist. For instance, for Q = 
circular orbits exist only for values of \L/{nM)\ > a/T2 ^ 3.45, whereas for Q = M and 




FIG. 2: The effective potential V for a neutral particle of mass ^ in a RN black hole of charge- 
to-mass ratio Q/M = 0.5, is plotted as function of r/M in the range [1.87,8], and the angular 
momentum L/fi in [0,10]. In this case the outer horizon r_|_ = 1.87M, and r^^ = 2.823M (see 
text). Circular orbits exist for r > 2.823M. The solid line represents the location of circular orbits 
(stable and unstable). The last circular orbit is represented by a point. The number close to the 
plotted point represent the radius r/M = 5.6, the angular momentum L/fi = 3.34 and the energy 
E/M = 0.939 of the last stable circular orbit. 
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FIG. 3: The effective potential V for a neutral particle of mass /x in a RN black hole of charge Q 
and mass M is plotted in terms of the radius r/M for different values of the angular momentum 
L* = L/(Mfj,) and Q = 0.5M. The outer horizon is located at r+ ~ 1.87M. The effective potential 
has a minimum, Vmin ~ 0.93, at rmin ~ 5.60M , for L* ^ 3.34. The plotted points represent local 
extrema. 

Q = 0.5M the existence condition implies that |L/(/xM)| > ^8 ^ 2.83 and |L/(/iM)| > 3.33, 
respectively. 

In this context, it is interesting to explore the stability properties of the circular motion 
at r = Tco- To find the explicit value of the last stable radius we solve the condition (|8]) in 
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FIG. 4: The radius of the last stable orbit tlsco/M is plotted as a function of the ratio Q/M. 
Numbers close to the points represent the value of the energy E/ ix and of the angular momentum 
L/{fiM) (underlined numbers) of the corresponding orbit. 



the black hole region and find 
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As expected, in the limiting case Q — )■ 0, we obtain the Schwarzschild value r^^co ~ The 
value of Tlsco decreases as Q/M increases, until it reaches its minimum value r'^^co ~ ^M at 
Q/M = 1. The general behavior of Tlsco in terms of the ratio Q/M is illustrated in Fig. HI 

Orbits with r > Tlsco are stable. Circular motion in the region r^^ < r < Tlsco is 
completely unstable. Since the velocity of a test particle at r = r^^ must equal the velocity 
of light, one can expect that a particle in the unstable region will reach very rapidly the 
orbit at r = Tlsco- For a static observer inside the unstable region, the hypersurface r = r^^ 
might appear as a source of "repulsive gravity" . This intuitive result can be corroborated by 
analyzing the behavior of energy and angular momentum of test particles. Indeed, Figs. [5] 
shows E/fi and L/{Mfi) as functions of r/M, for different values of the charge-to-mass ratio 
of the black hole. Both quantities diverge as the limiting radius r = r^^ is approached, 
indicating that an infinite amount of energy and angular momentum is necessary to reach 
r = r^^. As the ratio Q/M increases, the values of the energy and angular momentum at 
the last stable orbit decrease. For large values of the radius, the energy of circular orbits 
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FIG. 5: The pictures show (a) the energy E/ jj, and (b) the angular momentum Lj^/{M^) = L* of 
a circular orbit as a function of r/M. Notice that for Q = 0, the relevant radii are r_|_ = 2M and 
?^7+ = 3M, and the minima are Emin/fJ' ~ 0.943 and ~ 3.46, with r^m = 6M. Furthermore, 
for Q = 0.5M, we obtain r-|_ ~ 1.87M and r^^ ~ 2.83M, so that the minima are Emin/^^ ~ 0.939 
and L^j^^ ~ 3.34, with r^m — 5.61M. The energy and angular momentum diverge as the limiting 
radius is approached. 




FIG. 6: The energy E/ ^ of a test particle on a circular orbit as a function oi r/M and charge- 
to-mass ratio Q/M £ [0, 1]. The radius of the last stable circular orbit is also plotted (thick line). 
Numbers close to the points correspond to the energy and the angular momentum (underlined 
numbers) of the last stable circular orbit. 

approaches the limit E = fi, and the angular momentum L/{Mfi) increases monotonically. 
A more detailed illustration of this behavior in the case of the energy of the test particle is 
represented in Fig. |6] which shows E/ii in terms of the ratio Q/M and the radial distance 
r/M. 

The analytical expressions for the energy and angular momentum at the last stable orbit 
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FIG. 7: Plots of the energy Ei^sco/M (a), and the angular momentum Llsco/(m-^) (b) of the last stable 
circular orbit in terms of the ratio Q/M of the black hole. 

can be obtained by introducing the expression f ll2p into Eqs. (E]) and ([7]). A numerically 
analysis of the resulting expressions show the behavior depicted in Fig. ([7]). 

As a general result we obtain that the values for the radius of the last stable orbit as well 
as of the corresponding energy and angular momentum diminish due to the presence of the 
electric charge. Physically, this means that the additional gravitational field generated by 
the electric charge acts on neutral particles as an additional attractive force which reduces 
the radius of the last stable orbit. 




A. Extreme black hole 

In the case of an extreme black hole [Q = M) the outer and inner horizons coincide at 
r± = M. The effective potential vanishes at the horizon and tends to 1 as spatial infinity is 
approached. In this open interval no divergencies are observed. This behavior is illustrated 
in Fig. [8] where the effective potential 



is plotted for different values of the angular momentum L/{Mn). For Q = M the radius of 
circular orbits is 

The energy and angular momentum of test particles moving along circular orbits are given 
by 

_ {r-Mf L2 _ Mr2 

7^ ~ r2(r-2M)' 'J^ ~ r - 2M ' ^ ' 
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FIG. 8: The effective potential V for a neutral particle of mass /i in the field of an extreme RN black hole 
is plotted as function of the radius r/M for different values of the angular momentum L* = L/{Mii). The 
outer horizon is located at r+ = M . For L* k, 2.83 the effective potential has a minimum Vmin ~ 0.91 at 
fmin — 4M. There is a maximum Vmax ~ 1-35 at r « 2.19M for L* = 5, and for L* = 10 the maximum 
Vmax ~ 2.55 is located at r w 2.04M. 
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FIG. 9: Plots of (a) the energy E/fi and (b) the angular momentum L+/{Mjj,) = L* oi a. neutral test 
particle in circular motion around an extreme RN black hole. The outer horizon is r-|_ = M and r-,^ = 2M. 
At r = AM the minimum values of the energy Emin/ fJ^ ~ 0.91 and angular momentum ij^j„ ~ 2.83 are 
reached. 

Consequently, timelike circular orbits are restricted by the condition r > r^^ = 2M. Figure 
[9] shows the behavior of these quantities in terms of the radial distance. As r — )■ 2M, the 
energy and angular momentum diverge, indicating that the circular motion at r^^ is possible 
only along null geodesies. As expected, the local minimum of these graphics determine the 
radius of the last stable orbit. At r = rLsco = 4M, the energy is ~ 0.918/i and the angular 
momentum L = 2\/2nM. 
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(a) (b) 

FIG. 10: Plot of the energy E/ ^ of a neutral particle in circular motion around a RN naked singularity as 
a function of r/M and Q/M in the interval [1, 1.6] (a) and in the interval [1.1,2] (b). The line r* = Q^/M 
is also plotted. Numbers close to the points represent the energy and the angular momentum (underlined 
numbers) of the last stable circular orbits. 

III. NAKED SINGULARITY 

In the naked singularity case Q > M and the energy and angular momentum of the test 
particle can be written as 

E_ r^-2Mr + Q^ L+ _ I Mr - Q"^ 

~ - r^r^mr + 2Q^' TT ~ V r2 - 3Mr + 2Q2 " ^ > 

In Fig. [TD]a three-dimensional plot shows the energy as a function of both the circular orbits 
radius and the charge-to-mass ratio of the black hole. 

From the expressions for energy and angular momentum we see that it is necessary to 
consider four different cases: the value r = r^, = Q^/M, the region inside the interval 
1 < gVM^ < 9/8, the value QVM^ = 9/8, and finally the region defined by QVM^ > 9/8. 
In Figs. [miT2] the behavior of the effective potential is exemplified for different values of the 
ratio Q/M. 

Notice that for a RN naked singularity {Q/M > 1) the following inequality holds 

r* < < r^^ , (17) 

where r^^ = [3M ± ^J {9M'^ — 8Q'^)]/2 are the radii at which the value of the angular mo- 
mentum and the energy of the test particle diverge. 
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FIG. 11: The effective potential V for a neutral particle of mass /x in a RN naked singularity with Q/M = 1.1 

is plotted as a function of the radius r/M for different values of the angular momentum L* = L/{M ji). The 
classical radius r* = Q"^ /M = 1.21M is represented by a dashed line. For L* = the effective potential 
presents a minimum Vmin ~ 0.42 at rmin = ''*• For L* w 3.37 the minimum Vmin/ (J- ~ 0.95 is situated at 
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FIG. 12: The effective potential V for a neutral particle of mass in a RN naked singularity. Plot 
(a) is for Q/M = 1.5 so that r* = Q'^/M = 2.25M. For the lowest value L* = at r = r* there is 
global minimum with Vmin ~ 0.74. The value of Vmin increases as L* increases. Plot (b) corresponds to 
Q/M = 1.06, and shows the characteristic radii = 1.12M, r^^ = [3M + ^(QA'P - 8Q2)]/2 « 5.55M and 
= [3M - 1/(971^2 - 8Q2)]/2 « 1.44M. For L* = the global minimum is at r = n with Vmin « 0.33. 
As L* increases, the value of Vmin increases, and at L* w 2.7 a second local minimum appears. The first 
one with Vmin ~ 0.81 is at rmin ~ 1.29M and the second one with Vmin ~ 0.91 is located at rmin ~ 3.5M. 
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A. Static test particles 



Consider the orbit at r = r^, . In the naked singularity case the time-hke condition for the 
velocity Ug is satisfied only for r > r^, (see Eq. flA4p ). In the limiting case r = r^,, a time-like 
orbit with 

. , E I W , , 

-rr^Q^ '''' 

is allowed. It is interesting to note that r,, = Q'^ /M coincides with the value of the classical 
radius of an electric charge which is usually obtained by using a completely different ap- 
proach. This value appears here as the radius at which a particle can remain "static" with 
respect to an observer at infinity. This is an interesting situation which can be explained 
intuitively only by assuming the existence of a "repulsive" force. The above expression for 
the energy of the particle indicates that only in the case of a naked singularity a real value 
for the energy can be obtained. In fact, in the case of a black hole, the radius is situated 
inside the outer horizon so that r,, cannot reached by classical test particles. We conclude 
that the "repulsive" force can be the dominant gravitational force only in the case of a naked 
singularity. 

B. The interval 1 < Q'^/M'^ < 9/8 

In the first region, for 1 < Q'^/M'^ < 9/8 (see Fig. [T2b). time-like circular orbits can 
exist in the regions r* < r < and r > r^^. The boundaries r = r^^ correspond to 
null geodesies, as can be seen from the expression for the velocity along circular orbits as 
defined in Appendix |X] [ug^ = 1). This implies that there are two regions defined by r < r* 
and r G [r-y_,r^^] where no time-like particles can be found. This behavior is schematically 
illustrated in Fig. [131 In the limiting case Q^/M^ — )■ 1, the classical radius coincides with 
and therefore the only particle that can remain "static" on the classical radius must be 
a photon. 

C. The case Q^/M^ = 9/8 

For Q'^/M^ = 9/8, the exterior and interior photon orbits situated at and coincide. 
The effective potential behaves as illustrated in Fig. [Ml Local minima can be found in 
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FIG. 13: The pictures show (a) the energy E/^ and (b) the angular momentum L+/{Mii) = L* of a 
neutral particle moving on circular orbit around a RN naked singularity with Q/M — 1.06 < -^/O/S. The 
shaded regions arc forbidden, i.e., circular orbits can exist only in the regions defined by r^,, < r < r-y_ 



and r > r^^, where = Q'^/M w 1.12A/, = [3Af - ^ {'dAP - SQ'^)]/2 w 1.44M, and r^_^ = [MI + 



•y/ (9-M2 — 8(5^)]/2 « 1.55M (see text). A minimum for the energy and the angular momentum is located 
at r^m ~ 3.55Af, where L*^^^ 2.69 and Emin/fJ- ~ 0.91. Moreover, L*{r^,) — and [Emin/ y^ir*) ~ 0.33. 




FIG. 14: The effective potential V for a neutral particle of mass /i in a RN naked singularity with Q/A/ = 
•^9/8 is plotted as function of the radius r/M for different values of the angular momentum L* = L/^Mjj). 



In this case, = Q'^/M = 1.12A/ and r^^ = [3Af ± ^(9A//2 - 8g2)]/2 = l.SAf. 

different regions, depending on the value of the angular momentum of the test particle. 
Time-like circular orbits exists for all r > r,,, except at r = r^± = 3M/2, which corresponds 
to a photon orbit. The energy and angular momentum of circular orbits are given by 



E r2 - 2Mr + W 

o 



r-lM 



(19) 



11 r{r-lM^ At , 

and are plotted in Fig. [151 As the photon orbits are approached, the particle velocity tends 
to z/g = 1 and the energy and angular momentum diverge. 
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FIG. 15: Behavior of (a) the energy Ejy, and (b) the angular momentum L+/{Mii) = L* of a neutral 
particle moving along a circular orbit around a RN naked singularity with QjM = -^/O/S. Shaded regions 
are forbidden for time-like particles. 

D. The region Q'^/M'^ > 9/8 

For Q'^/M'^ > 9/8, the effective potential behaves as illustrated in Figs. [TT] and [T2k. 
Time-like circular orbits can exist for all r > r*. The time-like velocity condition {ug < 1) 
is always satisfied so that no limiting photon orbits can exist in this case. The energy and 
the angular momentum are plotted in Figs. [16] and [171 The angular momentum increases as 
the radius of the orbit r/M increases. In the limit of large values of r, the energy tends to 
E = fi. 

E. The last stable circular orbit 

To analyze the stability of circular orbits around a RN naked singularity we solve Eq. ([8]) 
under the assumption that Q > M. It turns out that real solutions exist only in the interval 
1 < Q/M < y/5/2. They can be represented as 



Isco 



Isco 



2M + 2 J4M2 - 3g2 cos 



1 /8M^ 
- arccos I 



2M-2J4M2 -3Q2sin 



arcsm 



and 



2M - 2 J4M2 - 3Q2 sin 



vr 1 

— I — arccos 
6 3 



M(4M2-3Q2)3/2 J 

M (4M2 - 3g2)3/2 j 
^8M^ - 9M2Q2 + 



M (4M2 - 3Q2)3/2 
However, it can be shown that in this interval it holds that Vc < Q'^/M 



(20) 
(21) 

(22) 
this 



solution is located inside the classical radius where no time- like circular geodesies are allowed. 
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FIG. 16: The pictures show (a) the energy E/ ^ and (b) the angular momentum L+/{M /i) = L* of a neutral 
particle moving along a circular orbit in a RN naked singularity with Q/M = 1.1 as a function of r /M . The 
shaded regions are forbidden. Circular orbits can exist only for r > sa 1.21Af . A minimum of the energy 
and the angular momentum is located at Tmin ~ 3.08M, where L*^^^ w 2.57 and Emin/lJ- ~ 0.90. Moreover, 
L*(n) = and {E^„Jfi)(r,) « 0.33. 
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FIG. 17: The pictures show (a) the energy E/^ and (b) the angular momentum L+/{Mii) of a test particle 
of mass /i in the field of a RN naked singularity as a function oi r/M and for different values of the ratio 
Q/M , satisfying the condition {Q/MY > 9/8 (see text). Circular orbits can exist only for r > r^, = Q^/M. 



Moreover, r^^^^ < r+^„ 



in the entire interval, except at Q/M = \/5/2 where r 



Isco 



Isco' 



For 



Q/M > V5/2 no solutions of Eq. ([8]) exist in the region defined by r > Q^/M so that for 



{Q/MY > 9/8 the last stable "circular" orbit is located precisely at r 



QVM. This 



situation is sketched in Fig. [18] where it can be seen that the energy of a particle located at 
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FIG. 18: The radius of the last stable circular orbit (solid line) of a neutral particle moving in the 
field of a RN naked singularity. The ratio Q/M varies in the interval [1,2.2] in plot (a), and in 
the interval [1, 1.15] in plot (b). The dashed curve represents the classical radius r* = Q"^ /M . The 
dotted curve corresponds to photon orbits with radius = [3M — ^ (9M'^ — 8(5^)]/2, whereas the 
dot-dashed curve denotes the photon orbits with radius r^^ = [3M + (9M^ — 8(5^)]/2. Shaded 
regions are forbidden. In the interval 1 < Q/M < 1.061 circular orbits can exist only for < 
r < (all stable) and r > r^_|_ (unstable for r^^ < r < r^g^^ and stable for r > r/\,^^). For 
Q/M > 1.061 and r > the region of stability divides into two separated regions: < r < r^^^^ 
and r > r'^^^^- The numbers close to the plotted points denote the energy E/jj. and the angular 
momentum L/{^M) (underlined numbers) of the last stable circular orbits. 

r = increases as the charge-to-mass ratio increases. In the interval 1 < Q/M < \/5/2, 
the energy and angular momentum of circular orbits situated on r = r^^^^ increase as the 
ratio Q/M increases. On the contrary, in the interval a/9/8 < Q/M < \/5/2, the energy 



and angular momentum of circular orbits with radius r 
decreases. This is due to the fact that at Q/M = 



Isco 



increase as the ratio Q/M 



9/8, the radius rf^^^ coincides with the 
radii r-y+ and which correspond to photon-like orbits. For Q/M > \/5/2, the last stable 
circular orbit is situated at the classical radius r^, = Q'^/M and the entire region r > r* is a 
region of stability. For ./9/8 < Q/M < V5/2, there are two regions of stable orbits, namely, 



< r < r 



Isco 



and r > r^"^^^, separated by a zone of instability defined by r^^^^ < r < r 



Isco' 



For 1 < Q/M < y9/8, two regions of stable orbits appear at < r < and at r > r^^^o! 
these regions are separated by the forbidden region, < r < r^^, in which no circular 



time-like orbits exist, and by the instability region located at r, 



7+ 



< r < r 



Isco' 



16 






(c) 

FIG. 19: The effective potential V for a neutral particle of mass in a RN naked singularity 
is plotted as a function of the radius r/M and of the angular momentum L* = L/{Mfi). Solid 
lines represent the radii of circular orbits. The last stable circular orbit radius is denoted by a 
point. Plot (a) is for Q = I.IM so that = 1.04M, r-y^ = 1.95M. The classical radius is 
= Q^/M = 1.21M, where L* = and Vmin ~ 0.42. Unstable circular orbits are in r < 3.94M 
and at r = 3.94M, the values are L* 2.81 and V/fi ^ 0.91. Plot (b) is for Q = 1.08M so that 
r^j_ do not exist and = Q'^/M = 1.16M with Vmini'''^) ^ 0.37. Stable orbits are in the interval 
r^, < r < 1.68M and at the boundary, r = 1.68M, the values are L* = 3.55 and V/n = 1.10. The 
interval 1.68M < r < 3.35Af corresponds to unstable orbits with the boundary values L* = 2.64 
and V/fi = 0.90 at r = 3.35M. Stable orbits are located at r > 3.35M. Plot (c) is for Q = 2M 
and = 4M. Stable circular orbits exist for r > 4M and at the boundary y(4M) = 0.86/i. 

The curves representing the radii of stable circular orbits are also showed in Fig. f lT9|) 
where the effective potential is plotted for selected values of the Q/M as a function of the 
angular momentum and the radius. 




17 



r/M 



10 



(a) (b) 

FIG. 20: The pictures show the energy E/^ (a) and the angular momentum L^/{^M) (b) of 
a test particle moving along a circular orbit around a RN source as a function of r/M and the 
charge-to-mass ratio Q/M in the range [0,2]. The curve r = r* = Q'^/M is also plotted. The 
numbers close to the points denote the energy and the angular momentum (underlined numbers) 
of the last stable circular orbits. 

IV. BLACK HOLE VERSUS NAKED SINGULARITY 

The main results obtained from the study of RN black holes and naked singularities are 
here summarized and compared in the plots of Fig. [20] - [22l 

In the three-dimensional plots of Fig. [20] the energy and the angular momentum ([7]) are 
shown as functions of the circular orbits radius and the charge-to-mass ratio Q/M, for both 
black hole and naked singularity cases. 

Fig. EH shows the circular orbits radius as a function of the angular momentum for 
different black holes and naked singularities. The cases of a Schwarzschild spacetime and 
the Newtonian limit are also plotted for comparison. Clearly, for large values of r/M all the 
curves converge to the Newtonian limit. Finally, the study of the circular orbits stability for 
neutral test particles in the RN spacetime is summarized in the Figl22] where the radius of 
the last stable circular orbits is plotted as function of Q/M. 

The location and structure of the stability regions for neutral particles around a black 
hole differ in a very strong way from the case of a naked singularity. Black hole sources 
are characterized by the existence of only one stability region which can extend from a 
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Q=M 



Q=1.061M - 



2 4 6 8 10 

L/(/i M) 

FIG. 21: The radius r/M of circular orbits for a neutral particle of mass /x in a RN geometry of 
charge Q and mass M is plotted as a function of the angular momentum L/{iJ.M) for different 
values of the charge-to-mass ratio Q/M m. the interval [0,2]. The radius of circular orbits for the 
Newtonian limit (dashed curve) and for the Schwarzschild spacetime (dot-dashed curve) arc also 
plotted. The plotted points represent the last stable circular orbits. Black curves correspond to 
stable circular orbits, and gray curves denote unstable circular orbits. The dotted line represents 
the classical radius r = r* = Q^/M. 

minimum radius risco £ [4M, 6M] up to infinity. In the case of a naked singularity there 
are two possible scenarios. If the mass-to-charge ratio satisfies the condition Q/M > V5/2, 
there exist only one stability region which extends from the classical radius r* = Q^/M up 
to infinity. The second scenario appears in naked singularities with charge-to-mass ratio 
within the range 1 < Q/M < \/5/2. In this case, we have two stability regions separated 
by either a region of instability or a forbidden region, where no time-like circular orbits are 
allowed, and a region of instability. 
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FIG. 22: The radius r/M of the last stable circular orbit of a neutral particle of mass /x in a RN 
geometry of charge Q and mass M is plotted as a function of the ratio Q/M in the interval [0, 1.5]. 
Here r^± = [3M ± ^(9M^ - 8(5^)]/2 and r± = M ± /M2Trg2 (ggg tg^t). In the Schwarzschild 
case, Q = 0, we find risco = 6M, r+ = 2M and r-y^ = 3M. For an extreme black hole, Q = M, we 
find 

I'lsco — 4M with r+ — ikf and r^_,_ — 2il<f. The numbers close to the plotted points denote the 
energy E/n and the angular momentum L/{iiM) (underlined numbers) of the last stable circular 
orbit. 

V. CONCLUSIONS 

In this work we discussed the motion of neutral test particles along circular orbits in 
the RN spacetime. We studied separately the case of black holes and naked singularities, 
emphasizing their differences. Our analysis is based on the study of the behavior of an 
effective potential that determines the position and stability properties of circular orbits. 
We also investigated in detail the behavior of the energy and angular momentum of the test 
particle in all possible configurations of black holes and naked singularities. 

We found that at the classical radius r = r^, = Q'^/M circular orbits exist with "zero" 
angular momentum. This means that a static observer situated at infinity would interpret 
this situation as a test particle that remains motionless (r = Q'^/M =const, =const) as 
time passes. This phenomena can take place only in the case of a naked singularity and is 
interpreted as a consequence of the "repulsive" force generated by the charge distribution. 

Black holes turn out to be characterized by a single zone of stability which extends from 
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a minimum radius risco up to infinity. The radius of the last stable circular orbit has its 
maximum value of risco — 6M in the Schwarzschild limiting case, and reaches its minimum 
value of risco — 4M in the case of an extreme black hole. The situation is completely 
different in the case of naked singularities. If the mass-to-charge ratio satisfies the condition 
Q/M > y/5/2, there exist only one stability region which extends from the classical radius 
r* = Q^/M up to infinity. In this case, the classical radius determines also the radius of the 
last stable orbit that could be, in principle, as large as allowed by the amount of charge Q 
that can be associated to a given mass M. If the charge-to-mass ratio is within the range 
1 < Q/M < -\/5/2, we found two stabihty regions that are separated either by a region of 
instability, in the interval a/9/8 < Q/M < V5/2, or a region of instability and a region 



non-connected regions in which test particles can remain orbiting forever. 
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Appendix A 

Let us introduce an orthonormal frame adapted to the static observers: 




AV2 



1 



d, 



(Al) 



r 



rsin^ 



with dual 



AV2 



dt, L)f 



r sin 6d(j) 



(A2) 



r 



Consider the tangent to a (time-hke) spatially circular orbit li" as 



u = r{dt + Cdct>) = 7 + 



where T and 7 are two normalization factors: 



^'-i-9tt-C9Hr' and ^ {1 - 



1 



21 



which assures that UaU"" = — 1; where ( is the angular velocity with respect to infinity and 
u is the "local proper linear velocity" as measured in the frame flAl|) . The angular velocity 
( is related to the local proper linear velocity by 



c 



so that r = 7/a/— Here ( and therefore also u are assumed to be constant along the u- 
orbit. As usually we limit our analysis to the equatorial plane; as a convention, the physical 
(orthonormal) component along -de, perpendicular to the equatorial plane will be referred 
to as along the positive z-axis and will be indicated by z and so = — e^. Here we focus 
attention to the time-like circular geodetic u± such that Vu±u± = 0, co-rotating (C+) and 
counter- rotating ((^_) with respect to the assumed positive (counter-clockwise) variation of 
the 0-angle respectively. It results: 



(Mr - Q2 



(A3) 



so that 



I Mr - g2 



79 



A 



1/2 



3Mr + 



(A4) 



l9 



\/-9tt 



(A5) 



with time-like condition Vg <1 satisfied. 

The particle's four-momentum is given by = yuu" where u*^ is given by flA5l) . The 
conserved quantities associated with temporal and azimuthal Killing vectors and are 
respectively 

V^9tt 

where E and L are the particle's energy and angular momentum respectively. Substituting 
the (lA4p we infer: 



E 



A 



r\ r 



M 



(A6) 



/i rVr" -3Mr + 2Q^' M/i ^/M ^/i^^^^3^^FTW 
these results are in agreement with the results obtained by the study of the effective potential. 
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